Plasmon-phonon Strongly-Coupled Mode in Epitaxial Graphene 
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We report the dispersion measurements, using angle-resolved reflection electron-energy-loss- 
spectroscopy (AREELS), on two-dimensional (2D) plasmons in single and multilayer graphene which 
couple strongly to surface optical phonon (FK phonon) modes of silicon carbide substrate. The cou- 
pled modes show discrete dispersion behaviors in the single and bilayer graphene. With increasing 
graphene layers on SiC(OOOl), a transition from plasmon-like dispersion to phonon-like dispersion is 
observed. For plasmon-like modes, the dispersion is strongly damped by electron-hole pair excita- 
tions at entering single-particle continuum, while phonon-like mode is undamped. In the region free 
of coupling, the graphene 2D plasmon exhibits acoustic behavior with linear dispersion with slope 
and damping determined by the Fermi surface topology. 

PACS numbers: 71.45.Gm, 73.20.Mf, 73.61.Wp 



The electronic properties of graphene, a two- 
dimensional (2D) monolayer of carbon atoms connected 
in a honeycomb lattice, has become the subject of much 
experimental and theoretical interest since it became 
practical to produce single sheets of this graphitic ma- 
terial Due to its high carrier mobility, graphene has 
the potential to substitute silicon in current silicon-based 
microelectronics, and the growth of large domain size was 
already reported on silicon carbide (SiC) wafer 0. The 
2D graphene sheet is insolated from insulating SiC sub- 
strate, sustaining a strictly 2D electron gas (2DEG) sys- 
tem where the collective oscillation (plasmon) can prop- 
agate along the surface. The wavevector-dependent dis- 
persion characterizes the dielectric response of 2DEG, 
which is important in every aspect of electronic appli- 
cation of graphene. Superior to optical excitation in 
momentum transfer, angle-resolved electron energy loss 
spectroscopy (AREELS) provides a unique tool to char- 
acterize the dispersion behavior at finite wavevector, 
where electron-electron correlation effects become evi- 
dent. The first observation of coupling between surface 
optical phonon and surface plasmon was on the surface 
of doped GaAs Q , the coupling strength was about 20 
meV. This dynamic property attracted much attention 
due to the wide application of GaAs in high electron mo- 
bility devices. 

We report experimental observation of the coupling be- 
tween the 2D "sheet plasmon" in graphene and "Fuchs- 
Kliewer" surface optical phonon (FK phonon) Q in SiC 
substrate. The "sheet plasmon" is collective oscillation 
of the 7r-valence charge density which propagates within 
the graphene sheet. This low-energy 2D-plasmon mode 
is different to the bulk 7r plasmon mode, which is a collec- 
tive oscillation of the tt- valence electrons in bulk graphite 
Q and coexisted in our spectrum measurements. The 
graphene 7r-plasmon mode is a strictly 2D mode whose 
energy decreases continuously to zero in the longwave- 
length (classical) limit q — ► 0. FK phonon arising from 
displacement of lattice ions is normally observed at polar 



semiconductor surfaces, such as GaAs and GaP Q. 

Our previous measurements 0] showed a discrepancy 
at small wavevector (q <0.05A _1 ); the graphene sheet 
plasmon energy deviating strongly from acoustic behav- 
ior. Similar measurements on the base 6H-SiC(0001) 
substrate showed the presence of a large-amplitude FK 
phonon mode in this same energy- wavevector range 
The suspicion was that the two modes might be strongly 
coupled. 

In an effort to resolve any coupling, we tuned the 
smallest possible incident electron energy which was high 
enough to excite plasmon and had a strong intensity for 
off-specular measurement, so as to increase the momen- 
tum resolution of our measurements. We observed an un- 
usually strong coupling about 130 meV between the sheet 
plasmon in the graphene and the FK phonon in the SiC 
substrate. We report these measurements here. Also, we 
compared the behavior with that of increasing number of 
graphene sheets, 1 to 5 monolayers (ML), grown on the 
same SiC(0001) substrate. The observed strong coupling 
carries over into each layer. However, subtle changes in 
the damping and dispersion are observed, due to a more 
complex Fermi surface, the topology and carrier density 
changing from strictly 2D graphene to 3D graphite. 

Graphene layers were prepared on a 6H-SiC(0001) 
crystalline wafer surface (0.1 f2 cm resistivity) by solid 
state graphitization The sheet structure has been 
well characterized and floating on a carbon-rich (6v3 x 
6^)^30° interfacial superstructure [HJ EH- AREELS 
measurements were performed at room temperature and 
in ultrahigh vacuum (2 x 10~ 10 Torr base pressure), on 
graphene layers of four different thicknesses and on a bare 
hydrogen-etched 6H-SiC(0001) surface. A medium en- 
ergy resolution of 10 meV was set to achieve strong signal, 
and the momentum resolution was better than 0.01 A -1 . 
[20I ] . which is small enough to resolve the concentrated 
peak profile of dipole scattering [7| . The wavevector q of 
each measurement at a particular scattering angle 9 S is 
calculated by following equation reflecting the conserva- 
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FIG. 1: (a) scattering geometry of AREELS experiments, 8i 
= 60° . (b) the extraction of loss peaks in spectrum taken at 
8 =1° on 1 ML graphene, which is comprised of two Lorentz 
peaks and background. Loss spectra of (c) 1.2 ML, (d) 2.3 
ML, (e) 3.3 ML and (f) 4.7 ML graphene dispersing with 
increasing scattering angle 8s- Black arrows indicate the po- 
sition of loss peaks after the subtraction of background and 
the deconvolution of peak envelope composing of two peaks 



tion of energy and momentum in plane, 



q = V2 ™° E [sinfl - y/l - E loss /E sin(0i - fl g )]. (1) 

where the loss energy E\ oss = E—E' and other definitions 
of variables are illustrated in Fig. QJa). The step between 
each measurement is determined by the impact energy E 
and the scattering angle 9 S which has a minimum step 
of 0.5°. A smaller step of wavevector can be achieved 



with smaller impact energy E, however the strength of 
spectrum signals also varies with impact energy E badly, 
so a considerable effort on tuning is needed to find the 
optimum setting to not only distinguish fine structures 
in q space, but also to yield enough signal strength at 
large off-specular scattering. In our work, the optimum 
setting for our spectrometer is E ~ 8 eV. 

Fig. Q] (c)-(f) show the energy loss peaks in 
the AREELS spectrum of low energy electrons back- 
scattered from graphene of various thicknesses on 
SiC(0001). The black arrows in Fig. [Hc)-(f) indicate 
the peak positions in loss spectra, which are the super- 
positions of multiple loss structures and the background, 
shown in Fig. LTJb). Generally the loss peaks disperse 
with increasing scattering angle 6 S , corresponding to in- 
creasing wavevector parallel to the surface. The surpris- 
ing spectroscopic structure is the splittings of peaks at 
certain energy, which is not disclosed in measurements 
with larger step at higher impact energy. The substrate 
of our samples, SiC, is a polarized ionic crystal with 
large band gap, and has attracted much attention in last 
decades, including researches on its surface excitation of 
various surface structures [12| • In order to clarify the ori- 
gins of this unexpected splitting of graphene dispersion, 
we took measurements on a bare 6H-SiC(0001) surface, 
which is different to 1 ML graphene only in top two layers, 
i.e. 6-\/3 buffer layer and top graphene layer, [ll[ so that, 
the polarization contributed from SiC substrate can be 
identified separately and used to analyzes measurements 
on graphene. 

Fig. HJa) shows AREELS spectrum taken on a bare 
6H-SiC(0001) surface before growing the graphene layer. 
Each spectrum of SiC show equally spaced peaks and the 
probability of multi-phonon excitation follows the Pois- 
son distribution, shown in Fig. (Hb), which is the charac- 
teristic of multi-phonon process [13[ . The calculated FK 
phonon frequency wfk is 116.7 meV in good agree- 

ment with experimental measurement of one-phonon ex- 
citation, 117 ± 10 meV and its multiples in the process of 
multi-phonon excitations. Under our experimental reso- 
lution (~ 10 meV), those FK-phonon peaks show little 
dispersion (straight dashed line in Fig. HI a)), in contrast 
to the large shifting of peaks on graphene layers (energy 
shifting A ~ 1 eV). As a comparison, Fig. [21a) also 
shows a normalized spectrum (black) taken on a single 
layer graphene in specular direction [6 S = 0°) at the same 
impact energy, where a drastic change of loss structure 
close to the elastic peak is due to the total reformation 
of atomic structure in the SiC(0001) surface. 

By using low-energy reflection electron energy loss 
spectroscopy with its shallow probing depth in the sur- 
face region, the collected loss structure in spectrum re- 
flects a long range dipole field decaying quickly in the 
bulk side. The excitation of FK phonons on SiC sur- 
face is actually the oscillation of induced field due to the 
displacement of positive and negative ions, which bears 
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FIG. 2: (a) loss spectra on SiC surface at various angles and 
normalized specular spectrum of 1 ML graphene (black). The 
loss spectra on SiC show multi-phonon loss peaks at 117 meV, 
236 meV, and 351 meV (magnified in insert), (b) the Poisson 
distribution of intensities of the specular spectrum 



a similar dielectric response properties to dipole field of 
2D plasmons induced by collective oscillation of electrons. 
This dielectric similarity enable the coupling in the region 
where the energy of FK phonons is comparative with that 
of 2D plasmon. But these two modes have very differ- 
ent origins and causes different screening behavior upon 
entering the single-particle excitation continuum, Fig. [3] 
Fig. [3] shows the dispersion behavior of the loss peaks 
observed in the AREELS spectra taken on graphene of 
various thicknesses. The distribution of single-particle 
excitations (SPE) due to intra- and inter-band transi- 
tion of 2DEG is also plotted with Fermi energy E-p = 
0.45 eV, and estimated Fermi wavevector fcp=0.079 A -1 
at electron density n = 2x 10 13 cm~ 2 [l(J. For 1.2 ML 
graphene (red dot in Fig. [3]), instead of smoothly in- 
creasing, the dispersion loses its continuity and jumps at 
q ~ 0.05 A -1 where crossovers occurs between plasmons 
and FK phonons (wfk >117 meV) . The insert in Fig. 
|3] shows a schematic change of dispersion curve where 
two modes couple with each other. Like the damping 
of plasmon dispersion of classical 2DEG, the cutoffs at 
entering SPE continuum are also observed in both sepa- 
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FIG. 3: The energy of Loss peaks plotted as a function of 
wavevector q calculated by Eqn. [T] The impact energy is 
around 8 eV. The dashed curves are guidelines for two dis- 
persion branches. Two shaded areas are single-particle exci- 
tations (SPE) continuum due to intra-band and inter-band 
transition. Insert: schematic of the coupling at which plas- 
mons disperse across FK phonon modes. 



rated branches, indicating that the coupled mode is dom- 
inated by plasmon feature, i.e. plasmon-like. Comparing 
with 1.2 ML graphene, 2.3 ML graphene (green square) 
exhibits similarly in splitting loss spectra and discrete 
dispersion. The slight difference comes from the second 
branch (0.05 < q < 0.1 A" 1 ), where the trend of data 
points is not as predictable as that of 1.2 ML, and can 
not be analyzed effectively if no more points are obtained 
with higher angular resolution of AREELS. The possible 
reason for this deviation is that electrons participating 
collective oscillation, which are also those electrons filled 
in Fermi level crossing tt bands, are affected by a splitting 
of 7r bands due to the formation of the second graphene 
layer 14j . 

With more graphene layers grown on SiC, the disper- 
sion exhibits dramatic change. In Fig. [4J the dispersion 
of 3.3 ML graphene(pink triangles) and 4.7 ML (brown 
diamond) graphene are plotted with that of single (red 
line) and bilayer (green line) graphene to show the dif- 
ference. Firstly, the energy and energy shifting of the 
coupled mode in 3.3 ML and 4.7 ML graphene arc much 
smaller than that of single and bilayer graphene. The dis- 
persions of 3.3 ML and 4.7 ML graphene are contained 
between one-phonon and two-phonon process and much 
flatter. Secondly, the dispersion continues into SPE, in 
contrast to the cutoffs observed in single and bilayer 
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FIG. 4: Dispersion behaviors of 3 ML (pink triangles) and 5 
ML (brown diamonds) graphene on SiC, in comparison with 
that of 1 ML (red line) and 2 ML (green line). 

graphene. This indicates a decreasing weight of plasmon 
feature in the coupled mode, consequently the disper- 
sion turns to phonon-likc which appears a much flatter 
dispersion and free of damping from electron-hole pair 
excitation in SPE. Accumulated with added graphene 
layer, the interaction between graphene layers increases, 
the characteristic of plasmon-phonon mode transits from 
plasmon-like to phonon-like. Basically this transition is 
due to the weakening of two-dimensional confinement of 
electron gas. 

Besides the discrete dispersion in \ow-q (q < 0.1 A -1 
) region due to plasmon-phonon coupling, another in- 
teresting observation is the linear dispersion in large-q 
(q > O.lA" 1 ). For an ideal 2DES, its plasmon frequency 
u>2D follows ^fq dispersion in the local-field approxima- 
tion in the longwavelength limit [ToT ]. Stern [id ] derived 
a higher-order correction including the effect of non-local 
field in a random-phase-approximation: 

, Anne 2 3 9 9 ,i 

- 2 D = [^ rT ^M + r | 9 2 + ..^ (2) 

where n, to*, and e s are the areal electron density, effec- 
tive mass and static dielectric constant of the medium, 
respectively. The first-order correction expresses in terms 
of the Fermi velocity v-p. At finite q, this correction term 
is dominant, plasmon disperse linearly. However, we note 
the recent discovery on metal surface, a so-called acoustic 
surface plasmon also shows a linear dispersion with [l7| : 

u>2D = av F q, (3) 



where coefficient a depends on the electron density of 
substrate and distance between 2D sheet and substrate. 
This linearity of plasmon is due to out-of-phase oscilla- 
tion between electrons in 2D sheet and in substrate. In 
our case, the epitaxial graphene bears a similar struc- 
ture, i.e. strictly 2D graphene sheet floating on a doped 
semiconductor substrate. We also notice the change of 
slope of dispersions for different graphene layers in large- 
q region, corresponding to the change of group veloc- 
ity of plasmon propagation, v-p. This is due to the 
scattering from degenerate ir bands induced by extra 
graphene layers. Furthermore, the change of slope is not 
monotonously decreasing, as the slope of 4.7 ML is larger 
than 3.3 ML in Fig. 01 indicative of the existence of bi- 
phase and out-of-phase polarization between layers which 
can be reflected in the oscillating charge distribution of 
multiple graphene layers [181 ] . The dispersion and damp- 
ing in this region (q > O.lA -1 ) also reflects the changing 
topology of the Fermi surface, evolving from 2D graphene 
to 3D graphite. 

In conclusion, we have observed an unusually strong 
coupling (~130 meV) between the dipolar electric fields 
generated by oscillations of the sheet charge of the 
7T electrons in graphene and the ionic charges in the 
substrate SiC. The sheet-charge density n ~ 10 13 cm~ 2 
is a result of charge transfer out of surface states on the 
SiC into empty 7r* states in the graphene sheet. Fuchs- 
Kliewer surface optical phonon modes in the SiC cause 
this charge transfer to oscillate. We believe this to be 
the origin of the unusually strong splitting observed. An 
usual feature, is that this coupling extends to graphene 
multilayers i.e. the multilayers are in noway screened 
by the first-layer, despite its very metallic character. 
This suggests charge transfer/scattering between layers. 
Further evidence is provided by changes in the plasmon 
dispersion of higher wavevectors q > O.lA -1 . However, 
the single-particle damping of the coupled mode does 
show different behavior, the dispersion continuing into 
the single-particle continuum in the case of the multilay- 
ers. This is indicative of subtle changes in momentum 
space, i.e. change in Fermi surface topology evolving 
from 2D graphene to 3D graphite, and changes in the 
screening properties. 
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